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The properties of low-energy neutron-proton systems are studied in an effective 
field theory where only nucleons figure as relevant degrees of freedom. With a finite 
momentum cut-off regularization scheme, we show that the large scattering lengths 
of the np systems do not spoil the convergence of the effective field theory, which 
turns out to be extremely successful in reproducing, with little cut-off dependence, 
the deuteron properties, the np 1 So scattering amplitude and most significantly, 
the MX transition amplitude entering into the radiative np capture process. 



1 Introduction 

In this talk I would like to report on a recent work on the application of effective 
field theory to two-nucleon systems carried out in collaboration with Kuniharu 
Kubodera, Dong-Pil Min and Mannque Rho. Effective field theories.^ (EFTs) 
have long proven to be a powerful tool in various areas of physicsjlra so it is 
no surprise that they are equally powerful also in nuclear physics. Indeed we 
have recently. harLlots of successful applications of EFTs in low-energy nuclear 
dynamics ffl'B'LfQ'ULl One of the spectacular cases was the chiral jperturbation 
theory (ChPT) calculation of the np — > dn process at thresholda with a per- 
fect agreement with experiment. In Ref. □, however, only the meson-exchange 
current correction relative to the one-body Ml transition amplitude was cal- 
culated, borrowing the latter from the phenomenological wave function of the 
accurate Argonne Vi$ potential tj. While thexe is nothing wrong there and in 
fact it is consistent with the strategy of ChPTH to calculate only the irreducible 
diagrams, there remained a "missing link" to a complete calculation, that is, 
to calculating everything within a given EFT. The motivation of this study is 
to do such a "first-principle" calculation M. In so doing we will compute the 
static properties of the bound np state (the deuteron) and the np scattering 
amplitude in the So channel. 

Limiting ourselves to the processes whose typical energy-momentum scale 
is much smaller than the pion mass, we keep only the nucleon matter field^ 
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c The anti-nucleon field is suppressed due to the largeness of the nucleon mass, and is also 
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as an explicit degree of freedom, integrating out all massive fields as well as 
the pion field. □ The resulting EFT is a non-relativistic quantum mechanics, 
with all the interactions appearing as a nucleon-nucleon potential in Lippman- 
Schwinger (LS) equation. The potential of the EFT is zero-range or contact 
interactions and their derivatives, because all the meson fields which mediate 
the nucleon-nucleon interactions are integrated out. This innocent looking sit- 
uation contains however many subtleties. First of all, even the leading order 
contact interactions are too singular to be solved by LS equation in three- 
dimensional space, thus we need to introduce a regulator and a renormalization 
scheme to handle the singularity. While the appearance of such a singularity 
is quite common in quantum field theories, the real subtlety comes from the 
fact that the np states in nature are all very close to the threshold: they are 
either weakly bound ( 3 Si) or almost bound ( 1 5'o)- Those states near thresh- 
old cannot be treated by perturbation expansion at all, which means that all 
the "reducible diagrams" up to infinite order should be summed by solving 
LS (or Schrodinger) equation. Furthermore, those states have huge scatter- 
ing lengths a, and the appearance of extremely small mass scales a -1 makes 
the convergence of EFTs by no means trivial. Indeed, using the dimensional 
regularization which has proven to be a very convenient and successful tool 
in handling singularities in most cases, Kaplan, Savage and WiseQ and Luke 
and ManoharQ have shown that the EFT breaks down at a very small scale, 

Pcrit = \J -^r for large scattering length a, where r e is the effective range and 

that this problem cannot be ameliorated by introducing the cionic degree of 
freedom. It was followed by the observation by Beane et alB that, for non- 
pcrturbative cases such as ours, the physical results after the renormalization 
procedure roay still depepd on the regularization scheme. As pointed out by 
Beane et alu and Lepagel^3 the problem can however be resolved if one uses a 
cut-off regularization. In EFTs, the cut-off has a physical meaning and hence 
should be set by the mass of the lightest degree of freedom which is integrated 
out, namely the pion in our case. If one chooses too low a cut-off, the valid 
region of EFTs unnecessarily shrinks down, while if one chooses too high a 
cut-off, one introduces irrelevant degrees of freedom and hence makes the the- 
ory unnecessarily complicated. We find that the optimal cut-off in our case is 
A ~ 200 MeV as one can see from the results in Table || and Figures I and 2. 



integrated out as is in heavy-baryon formalism. 
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2 Renormalization and Phase shift 

To do the calculation algebraically, we choose the following form of regulariza- 
tion appropriate to a separable potential given by the local Lagrangian: 

(p'|^|p} = Sa(p'V(p'-p)Sa(p 2 ) (1) 

where <Sa(p 2 ) is a regulator which suppresses the contributions from |p| > A, 
limipi^A Sa(p 2 ) = 1 and um | P |»A Sa(p 2 ) = 0, and V(q) is a finite-order 
polynomial in q. We shall do the calculation to the next-to-leading order 
(NLO), and so the most general form of V^(q) is 

Air 

V(q) =^{Co + iG^ + D 2 cxV)q*q>) , (2) 

where M is the nucleon mass and cr 1 - 7 is the rank-two tensor that is effective 
only in the spin-triplet channel, 

i „3 i „3 „i xii \ 

(3) 

Note that the coefficients Co, 2 are (spin) channel-dependent, and that D 2 is 
effective only in spin-triplet channel. Thus we have five parameters; two in 
1 So and three in 3 Si channel, which will be fixed from experiments. Since the 
explicit form of the regulator should not matterli 2 ] we shall choose the Gaussian 
form, 




S A (p^)=exp^-^j (4) 

where A is the cut-off. The LS equation for the wavefunction \tp), = 
\ip)+G° V\tp) where \ip) is the free wavefunction and G° is the free two-nucleon 
propagator depending on the total energy E, (p'\G°\p) — — ^ P J P ^ — leads to 

E— — +i0+ 

the S-wave function (for the potential (|2|)) of the form 

S A (ME)C E 
V'(r) = <p(r) + r „ 

1—1 E^E 



1-^{V 2 + ME) 



ZD 2 S 12 (f) d 1 d 
r- 



C E VS drr dr 
where Si 2 (f) = 3<Ji ■ r a 2 • f — o\ ■ a 2 , 



f A (r) (5) 
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with 



C E = a A U + ^a A r A ME] + (yZD 2 ME^Y E , 



a A = Z [Co + (Cl + 5 s .iDl)h)] , 
2Z 

r A = — [2C 2 - (Cf - Ss,iDl)h] 



(9) 

(10) 
(11) 



where I n (n — 2, 4) are defined by 



In 



dxx n Sl(x 2 h 2 ) 



(12) 



With the regulator (||), the integrals come out to be h = ~^/^^ 3 an d li = 

The phase shifts can be calculated by looking at the large- r behavior of the 
wavefunction. To do this, it is convenient to separate the pole contributions 
of the integrals Eqs.((|, 0) as 



T E = -tVMESl(ME) + I A (E), 



f A (r) 



Sa(ME) 



AVMEr 



H(Ar, 



ME, 
"A 5 "' 



(13) 
(14) 



which define the functions I\{E) = AI(^-) and H(kr, ^-), both of which are 
real. Note that H(Q, e) — 1 which makes Ta(0) finite, and that lirri x >i H (x, e) = 
0. The 1 Sq phase shift (5( 1 S'o) takes the form 



p cot 6C So) = 



1 



S\(ME) 



Ia(E) - 



1 



a A (l + ia A r A ME) 



(15) 



The "effective" low-energy constants, a A and r A , are fixed by comparing (15) 
to the effective-range expansion 



1 1 2 
pcot o = h —r e p + 

a 2 



(16) 
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A/(0) = 

2J'(0) 
A 



1 


A 


a 
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a 


dp 2 



p 2 =0 



2 

aA?" 



(17) 
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They then give us the "renormalization conditions" of the Co and C%, with a 
given value of A. Two important observations to make here: (a) We note that 
there is an upper bound of A, Am^x, if one requires that Z be positive and 
that C2 be real. That is, for A > AM ax , the potential of the EFT becomes 
non-Hermitian. With a — —23.732 fm_and r e = 2.697 fm for the 1 Sq channel 
taken from the Argonne w 18 potential^ (which we take to be "experimental"), 
we find that AMax — 348.0 MeV; (b) the value Az=i defined such that Z = 1 
when A = Az=i — 172.2 MeV is quite special. At this point, we have r A = 
and C2 = 0, that is, the NLO contribution is identically zero. This corresponds 
to the leading-order calculation with the A chosen to fit the experimental value 
of the effective range r e . A similar observation was made by Beane et alB using 
a square- well potential in coordinate space with a radius i?, with playing 
the role of A. 

The resulting phase shift with A = Az=i is plotted in Fig. 1. We see 
that the agreement with the result taken from the Argonne v\% potential E3 
is perfect up to p ~ /2. Beyond that, we should expect corrections from 
the next-to-next-order and higher-order terms. In Fig. 2, we show how the 
phase-shift for a fixed center-of-mass momentum, p = 68.5 MeV varies as the 
cut-off is changed. The solid curve is our NLO result, the dotted one the LO 
result (with C 2 =0), and the horizontal dashed line the result taken from the 
vis potential ("experimental"). We find that our NLO result is remarkably 
insensitive to the value of A for A > m^. It demonstrates that, going to the 
higher-order calculations in EFTs, we have not only more accurate agreement 
with the data but also less dependence on the cut-off. 

As for the 3 Si coupled channel, the phase shift 3 5i and the mixing angle 
ei are given as 
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FAME) 



Ia(E) - 



p cot 5( 3 5i 

V(E) 

1-V 2 (E) aA{l + \a K TKMEy 



l-r,\E) 



a A (l + ±a A r A ME) 



ZD 2 ME 



(19) 
(20) 



where T)(E) = — tanei and we have used the eigenphase paxametrizationlli 
The D 2 can be fixed by the deuteron D/S ratio rj d ~ 0.025 B at E = -B d 
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p (MeV) 

Figure 1: np 1 So phase shift (degrees) vs. the center-of-mass (CM) momentum p. 
Our theory with A = Az=i — |172 MeV is given by the solid line, and the results 
from the Argonne vig potential Ej ("experiments") by the solid dots. 
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Figure 2: np 1 So phase shift (degrees) vs. the cut-off A for a fixed CM momentum 
p — 68.5 MeV. The solid curve represents the NLO result, the dottedreurve the LO 
result and the horizontal dashed line the result from the wis potential cj. 
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Table 1: Co, C2 and D2 for various values of A. The unit of Co is GeV 1 while that of C2 
and D2 is GeV -3 . D2 in Sq channel are identically zero. 





1 Sq 


3 Si 


A (MeV) 


Co 


c 2 


C 


c 2 


D 2 


150.0 


-9.877 


-56.4 


-15.503 


-206.6 


186.9 


172.2 


-9.482 





-12.051 


-78.0 


169.7 


200.0 


-10.185 


37.9 


-9.925 


2.5 


156.3 


216.1 


-11.377 


51.8 


-9.730 


31.7 


152.7 


250.0 


-17.167 


74.2 


-14.342 


83.4 


158.1 



with Bd the binding energy of the deuteron, 



ZD, 



«A 



1 - —a\r\M Bd 



(21) 



The renormalization procedure is the same as for the 1 Sq channel. The only 
difference is that the value of Az=i that makes Z = 1 does not coincide-svith 
A rA= o that makes r\ — 0. Using a — 5.419 fm and r e = 1.753 fmliJ for 
the 3 Si channel, we find that A M ax = 304.0 MeV, A z =i = 198.8 MeV and 
A rA= o = 216.1 MeV. The values of Co, C 2 and D 2 with respect to various 
values of A are listed in Table [j] One can see the "naturalness" of the coefficient 
by making them dimcnsionless, for example, w^Co, rn^C 2 and m^D 2 in spin- 
triplet channel with A = 216.1 MeV are —1.4, 0.1 and 0.4, respectively. 



3 Results and Discussion 

Given Co, C 2 and D 2 for a given A, all other quantities are predictions. The 
binding energy of the deuteron is determined by the pole position, 



7 5i(- 7 2 ) + / A (- 7 2 ) 



a A (l - ia A r- A 7 2 ) 



(22) 



with 7 = \J MBd- The S- and D-w&ve radial wavefunctions of the deuteron 
are 



u(r) = c 
w(r) = r\d 



7r - H(Ar, 

r 2 d 1 d 1 
7 2 dr r dr r 



-7 
A 2 



a 47rr ^(3)/ \ 



e - 7r - H(Ar, 



-7 
A 2 



2 1 



(23) 
(24) 
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where 



(V^-l)A 2 m) 
fA a A (l + ±a A r A ME)I 2 S A (ME)- 1 1 

We now have all the machinery to calculate the deuteron properties: the 
wavefunction normalization factor ^4 S , the radius r^, the quadrupole moment 
Qd and the .D-state probability Pp, which are defined as 

/>oo 

A- 2 = / dr [u 2 (r)+u; 2 (r)] , 



(i 



4 

' |u(r)w(r) - — jJ- 



dr r 



P D =A Z S / druo 2 {r). (27) 
Jo 

The magnetic moment of the deuteron /i^ is related to the Pjj through 

= Ms ~ I (lis - I) Pd (28) 



2 V 2, 

where fis — 0.8798 is the isoscalar nucleon magnetic moment. Finally the 
one-body Jsovector Ml transition amplitude relevant for n + p — > d + 7 at 
thresholdca is 

/•OO 

M 1B = / dru(r)u {r) (29) 
Jo 

where uo(r) is the 1 S'o radial function, 



«o(r) = . , ag ; H{Ar,—-)+p A —Sy{r) 30 

smd( i ,bo) A z 

The (parameter-free) numerical results are listed in Table || for various 
values of the cut-off A. We see that the agreement with the experiments 
(particularly for A = 216.1 MeV) is excellent with very little dependence on 
the precise value of A. It may be coincidental but highly remarkable that even 
the quadrupole moment, which (as the authors of Ref. 1121 stressed) the vi$ 
potential fails to reproduce, comes out correctly. 
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Table 2: Deuteron properties and the Ml transition amplitude entering into the np capture 
for various values of A. 



A (MeV) 


150 


198.8 


216.1 


250 


B d (MeV) 


1.799 


2.114 


2.211 


2.389 


A s (fm-*) 


0.869 


0.877 


0.878 


0.878 


r d (fm) 


1.951 


1.960 


1.963 


1.969 


Qd (fm 2 ) 


0.231 


0.277 


0.288 


0.305 


Pd (%) 


2.11 


4.61 


5.89 


9.09 


Ud 


0.868 


0.854 


0.846 


0.828 


Mib (fm) 


4.06 


4.01 


3.99 


3.96 



Exp. 



2.225 



0. 

1.966(7) 
0.286 

0.8574 



IT 



2.225 
0.885 
1.967 
0.270 
5.76 
0.847 
3.98 



Let us compare our result ( |l5| ) with that obtained with the dimensional 
regularizationEl, 

P cotS \ Dlm . = 77— l 77™- ( 31 ) 

a(l + ^ar e Mhi) 



Expanding p cot S of (31 ) in ME, we find that the coefficient of the n-th order 
term is order of a™ -1 ?'". This increases rapidly with n when a is large, dis- 
agreeing strongly with the fact that the low-energy scattering is well described 
by just two terms o£-the effective range expansion in (filf). This observation led 
the authors of Ref. □ to conclude that the critical momentum scale at which 
the EFT expansion breaks down is very small for a very large a: 

Pcrit\ Dim ~ \ -^~- (32) 

V ill e 

We arrive at a different conclusion. With the cut-off regularization, the scat- 
tering length a is replaced by an effective one, aA,,-that is order_of A^ 1 for 
large \a\. This agrees with the findings of Beane et alia and LepageE3 Counting 
r e to be order of A -1 , the n-th order coefficient now is A 1_2n , as one would 
expect on a general ground. Recalling that a small (and negative) scattering 
length corresponds to weak (and attractive) interactions, it is also remarkable 
that the mechanism of the finite cut-off EFTs is quite similar to the "quasi- 
particle" phenomena since both convert highly non-linear systems into weakly 
interacting systems. 

Using the finite cut-off regularization scheme which is the most faithful 
way to realize the principles of EFTs, we have demonstrated that the low- 
energy nuclear physics can be well-described by EFTs. In particular, it is 
satisfying that the classic np capture process can be completely understood 
from a "first-principle" approach. Here the cut-off regularization was found 
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to be highly efficient. With the dimensional rcgularization the Ml matrix 
element was found to be in total disagreement with the result of the Argonnc 
vis potential. 

There are several important and urgent extensions and applications of the 
finite EFTs in nuclear physics. One of the most important applications is the 
proton fusion process, p + B-^ d + e + + i>, at threshold which plays a crucial 
role in the stellar_evolutionEJ This process has been recently worked out and 
will appear soonEj From the theoretical side, the most urgent task is to take 
into account the pion degree of freedom so as to extend the calculation to 
higher chiral order. This would enable us to study the interplay between the 
breakdown of an EFT and the emergence of a "new physics." 
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